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Abstract: In this paper we conduct a numerical study of the supersymmetric 0(3) non- 
linear sigma model. The lattice formulation we employ was derived in Q and corresponds 
to a discretization of a twisted form of the continuum action. The twisting process exposes 
a nilpotent supercharge Q and allows the action to be rewritten in Q-exact form. These 
properties may be maintained on the lattice. We show how to deform the theory by the 
addition of potential terms which preserve the supersymmetry. A Wilson mass operator may 
be introduced in this way with a minimal breaking of supersymmetry. We additionally show 
how to rewrite the theory in the language of Kahler-Dirac fields and explain why this avenue 
does not provide a good route to discretization. Our numerical results provide strong evidence 
for a restoration of full supersymmetry in the continuum limit without fine tuning. We also 
observe a non- vanishing chiral condensate as expected from continuum instanton calculations. 
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1. Introduction 

Supersymmetric field theories are an integral part of modern theories of particle physics. 
They provide a framework for solving the gauge hierarchy problem [g] by eliminating many 
divergences typical of quantum field theories through the cancellation between fermionic and 
bosonic loops. Moreover in the large N limit, it is known ||, Q that supersymmetric gauge 
theories are related to quantum gravity and string theory. Two dimensional sigma models on 
the other hand are important because they have a rich mathematical structure, and moreover 
there exists a deep rooted analogy between them and four dimensional Yang Mills theories |^]. 
Indeed the former can serve as perfect theoretical laboratories to test methods and approaches 
developed for solving the problems of these far richer and more complicated theories. For 
example, the magnetic monopole and dyon in the gauge theory correspond to the kink and Q- 
kink solution of sigma models. In many cases the interesting physics lies in non-perturbative 
regimes which motivates use of the lattice to study these systems. 
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Unfortunately generic discretizations of super symmetric theories break supersymmetry 
explicitly and necessitate fine tuning the couplings to a (usually large) number of lattice oper- 
ators in order to approach a supersymmetric continuum limit. Recently, attention has turned 
to formulations which aim to preserve one or more exact supersymmetries at non-zero lattice 
spacing 1 . The hope is, this residual supersymmetry will shield us from the appearance of 
relevant operators in the lattice effective action which violate the full symmetry group thus 
reducing or possibly even eliminating such fine tuning problems. Two distinct approaches 
in this direction have been followed; in the first pioneered by Kaplan and collaborators, the 
lattice theory is constructed by orbifolding a certain supersymmetric matrix model and then 
extracting the lattice theory by expansion around some vacuum state @, |, |, [fij. This has 



recently been extended to the the interesting case of gauge theories in two dimensions with 
matter fields interacting via a superpotential in [|llj . The second approach relies on discretiz- 
ing a twisted version of the supersymmetric theory. Twisting was first introduced by Witten 
] |l^ ] in the context of topological field theories. It consists of constructing a new rotation 
group from a combination of the original rotation group and part of the i?-symmetry associ- 
ated with the extended SUSY. The supersymmetric field theory is then reformulated in terms 
of fields which transform as integer representations of this new rotation group [13, 14, 15, 16]. 



In flat space one can think of the twisting as a merely exotic change of variables in the the- 
ory. When applied to the supersymmetry algebra, a scalar nilpotent supercharge is exposed. 



Furthermore as argued in [17, |18| the twisted superalgebra implies that the action rewritten 
in terms of these twisted fields is generically Q-exact. In this case it is straightforward to con- 
struct a lattice action which is Q-invariant provided only that we preserve the nilpotency of Q 
under discretization. An other approach in the same direction ][l5|, 19], attempts to preserve 
all the twisted supercharges on the lattice by introducing a non-trivial commutation relation 
between the coordinates in superspace. However, this formulation remains controversial after 
a recent paper |2(| pointed out an inconsistency appearing in this method when applied to a 
toy quantum mechanical model. Finally, although the orbifolding and the twisting approaches 
appear different, Unsal pi]] recently showed that the former approach reproduces the twisted 
Yang-Mills theories in the continuum limit. 

The twisting approach to constructing lattice theories was initially developed for theo- 
ries without gauge symmetry |2^, |2j|. An implementation for supersymmetric lattice gauge 
theories based on balanced topological field theories was given by Sugino An approach 
emphasizing the geometrical nature of the twist and employing Kahler-Dirac fermions was 



then developed by Catterall |17|, 25] to construct super Yang-Mills theories in two and four 
dimensions. The use of Kahler-Dirac fermions for formulating lattice supersymmetry was first 
proposed in [p6|. 



2. The 2D continuum action 

As was shown in |ffl|, the action of a general two dimensional sigma model with Af = 2 



1 We primarily discuss Euclidean formulations - for recent work on Hamiltonian approaches see |s| 
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super symmetry may be written, using complex coordinates, in the twisted form 



S = (3 j d 2 a (2h + -g I jd + <j) I d^ J 

- h+'g^D^ 7 - ^-gjjvlD^ + h + ~ R nn vW-^ J V> 7 ) (2.1) 



with 



D + ^ J = d^ J + T J KL d+4> K i> L (2.2) 



The requirement of M = 2 supersymmetry forces the target manifold to be Kahler in which 
case gu = gjj = and the only non-zero Christoffel symbols are Tj K and Ij^. (see Appendix 
[A]) . This action is invariant under four supersymmetries as expected for a theory with Af = 2 
supersymmetry in two dimensions. In the twisted construction we focus on a single hermitian 
twisted supercharge Q whose action on the fields is given by 



(2.3) 



= if; 1 Qcj) 1 = if; 1 

Q^j 1 = Qip 1 = 

Qr)\ = B{- T^j^rif Qr]L = bL- t rj ^ 3 ^ 

QB{ = -T^B* - R I KjL il> K il> J rfi QBi = -T T JR ^ J B? - R 1 - JLR ^^ h r^ 

The action of the remaining charges of the N = 2 twisted supersymmetry is given in appendix 
0. The field B\ is an auxiliary field introduced to linearize the transformations and render 



the transformation nilpotent off-shell. It has been removed from the action in eqn. (2.1) 
by employing its equation of motion B^_ = d+cj) 1 . The invariance of the action under this 
supercharge Q follows just from the nilpotency of the latter and the fact that the above action 
is Q-exact - that is it can be written as the Q-variation of some function which, borrowing 
from BRST gauge fixing terminology, can be termed a gauge fixing fermion O. 



3. The 2D lattice action 

Surprisingly translating the action in ([O]) to the lattice is pretty straightforward. Indeed, as 
the twisted Q-symmetry makes no reference to derivatives of the fields its nilpotent property 
is preserved when continuum fields indexed by a continuous coordinate are replaced by lattice 
fields carrying a discrete index. We then obtain the supersymmetric lattice action by just 
replacing the continuum derivative by a symmetric finite difference 

s = /?^(2/ 1 +- 5/7 a;/ai<a 7 

X 

~ h+'g^D^ 1 - h + -gj jV 7 _D + i; J + ^h+~ R i7 jjvW-^ J V 7 ) (3.1) 

where now, 

D + ^ J = A s + ip J + Y J KL A\$ K i) L (3.2) 
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and Af 



— A* rb zA| where A* = ^(A+ + A~) and A 1 * 1 are the usual forward and backward 
difference operators. However the lattice action in the form shown in eqn. (|3.1| ) suffers from 
the usual fermion doubling problem. Indeed it is easy to see that the kernel of the (free) 
lattice Dirac operator constructed this way contains extra states which have no continuum 
interpretation. For instance consider the fermionic part in (|3.1|). If we ignore the interaction 
term in (B.2I) then one finds that in Fourier space the kernel of Ai corresponds to solving 







e tkx [sin ki + i sin fe 2 ] 



(3.3) 



Besides the solution (fcijfe) = (0,0), one also has (0, 7r) , (it, 0) and (tt,7t). The doubling 
problem is made worse because supersymmetry propagates it to the bosonic sector. There 
are two obvious lines of approach one may take to avoid this problem. One is rewrite the 
continuum theory in the language of Kahler-Dirac fields and try to utilize the discretization 
procedure developed in [[l?]] to construct the lattice theory. This approach is reviewed in the 
next section where it is shown that such a procedure runs into difficulties. 

The obvious alternative, which we have pursued here, is to add some form of Wilson 
mass term to lift the mass of the doubles up to the scale of the cut-off. However an ad hoc 
addition of such a Wilson term will break supersymmetry explicitly in an uncontrolled way 
thus spoiling our goal of preserving an element of SUSY on the lattice. Fortunately, it was 
shown |2{| that in the case where the Kahler manifold possesses some isometries, it is 
possible to add potential terms that involve the holomorphic Killing vectors V 1 associated 
to those isometries in such a way as to keep the action Q-exact. By a careful choice of such 
potential terms we may add Wilson operators which accomplish the task of rendering the 
doublers heavy. In complex coordinates the possible terms are, 



AV j V/ + X^'VjVji/j" - ^h + -r ] 1 _VjVjri J + + h.c 



1 



— »J ' 



(3.4) 



Here, A 2 is an arbitrary parameter. Here, we choose A = h + = ^. To keep the action 
Q-exact the action of the twisted supersymmetry must be changed. 



Q^ 1 = ip 1 
Q^ 1 = V 1 



(3.5) 



Q is no longer nilpotent but Q 2 just amounts to a Lie derivative with respect to the Killing 
vector field. Indeed one can easily show that, 



Q 2 Bi 



v 1 

djV 1 ^- 7 

djv 1 ^ 

djV l B J + 



(3.6) 
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To introduce a mass term would require utilizing a Killing vector of the form 



V 1 = im^ 1 (3.7) 

It should be clear that the transformation induced by such a Killing vector corresponds to a 
infinitessimal phase rotation on the complex fields and can be seen to be a good symmetry 
of both the continuum and lattice actions. Many Kahler manifolds possess metrics invariant 
under such a phase rotation, for example, the CP N theories and specifically the CP 1 ~ 
0(3) model considered in detail later. In the latter case the theory actually possesses three 
isometries corresponding to the global 0(3) symmetry of the theory - the phase rotation 



associated with the Killing vector eqn. 3/7 corresponds to invariance under rotations about 
the z-axis. 

To remove the doubles in lattice regularizations of such models we have employed such a 
twisted mass term with the constant mass parameter replaced by a Wilson operator m — > mw 
where 

m w = \ (A+AI + A±A;) (3.8) 
This leads to the additional terms in the action 

AS = ^2 ^mwtfmwfa + ^ 2 i^Vjg J -jrnw<t> J '4> J — ^h + ~rjLX7jg J jmw4> J 'n'l + h- c 

(3.9) 

The generalized phase rotation associated with mw n o longer corresponds to an exact isom- 
etry of the metric and hence the modified action, while still Q-exact, is no longer exactly 
invariant under the generalized Q-symmetry. Nevertheless we will argue that the twisted 
Wilson operator acts as a soft breaking term and hence should not affect the renormalization 
of the lattice theory for small enough lattice spacing. 

The argument goes as follows. In the continuum one is used to thinking of mass terms 
as serving only to break supersymmetry softly. This expectation relies on the idea that any 
mass parameter will be small compared to the U.V cut-off in the theory. However, the generic 
Wilson operator in a lattice theory does not satisfy this property since the potential doublers 
in the theory pick up a mass from the Wilson term on the order of the cut-off. Thus addition 
of such Wilson terms generically will lead to a hard breaking of supersymmetry and it will 
be necessary to fine tune additional operators to recover a super symmetric continuum limit. 
However, the twisted mass operator we use here does not have this property - in the limit 
of small a the propagators are dominated by contributions from a state near the origin of 
the Brillioun zone and a set of would be doublers. The latter contribution does not break 
supersymmetry since the doublers contribute like additional states with a constant twisted 
mass. Such a mass deformation preserves the Q-symmetry of the lattice action. 

The soft character of the breaking can be seen in yet another way. For small lattice 
spacing a (large (3) we can expand the bosonic action to quadratic order. Subsequently 
integrating over the bosons yields a determinant dei(m w — D^_D^_). But at one loop this 
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determinant is cancelled by an identical fermionic contribution det(M) where 



M 



-D* 



-D_ mw 



Thus the lattice theory appears both double free and supersymmetric at large (3. These 
arguments lead us to expect the Q-symmetry to be restored at large coupling (3 without 
additional fine tuning. If this is the case then simple power counting arguments lead us to 
conclude that the only (marginally) relevant counter terms must take the Q-exact form 



= Q{v ia fij(<t>) {d a ^+B j )) 



(3.10) 



where we have reverted to a formulation using real fields. General covariance ensures that 
fij is a tensor which may then be taken to represent a quantum renormalization of the 
target space metric tensor. This counterterm structure would be consistent with a lattice 
model which exhibits M = 1 supersymmetry in the continuum limit. The restoration of full 
J\f = 2 supersymmetry appears to require additional constraints. Luckily, such constraints 
are present in the form of additional discrete symmetries of the lattice action. Consider the 



classical action in Kahler form given in eqn. 3.1. It is trivial to see that this action is also 
invariant under the finite transformations 



—irjL 



(3.11) 



Actually, this additional symmetry arises from the Kahler structure of the target space ap- 
pearing in the classical action j30|. This additional symmetry of the lattice model then ensures 
that only counterterms compatible with a Kahler target space survive in the quantum effective 
action. But as was shown in any model with M = 1 supersymmetry and a Kahler target 
space automatically possesses M = 2 supersymmetry. Thus on the basis of these arguments 
we expect that no additional fine tuning is needed to regain the full supersymmetry of the 
continuum model which as will be seen later, is confirmed by our numerical results. 



4. Kahler-Dirac reformulation 

Recent work has emphasized the geometrical nature of certain twisted super Yang-Mills theo- 



ries 



17, 25 1 . In these constructions all fields appear as antisymmetric tensors with the twisted 



fermions arising as components of a geometrical object called a Kahler-Dirac field. The ar- 
guments leading to this geometrical interpretation are quite general, depending only in the 
number of supercharges and i?-symmetry, and imply that it should be possible to formulate 
the twisted sigma models also in this language. Furthermore, the geometrical construction 
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is key to a successful discretization of the Yang-Mills theory avoiding fermion doublers [27]. 
Thus such a reformulation of the sigma model potentially offers another way to build a su- 
persymmetric lattice theory without encountering fermion doubling. As such it would offer 
an alternative to the addition of Wilson operators as described in the previous section. As we 
will see this reformulation can indeed be done in the continuum but an obstruction prevents 
any simple translation of the continuum theory to the lattice. 

As shown in, for example, p5| the basic M = 2 twisted supermultiplet in two dimensions 
involves two Kahler-Dirac fields, a bosonic one = ((f), An, -B12) whose p-form components 
commute and fermionic one = (A, 77^X12) with grassmann valued forms. Initially we 
will consider a flat target space and consider only the continuum theory. The action of the 
nilpotent twisted supersymmetry is simply 



= A QX = 
Qrjf, = QA^ = 

QB 12 = Qxn = 



(4.1) 



These transformations are essentially the same as the corresponding Yang-Mills variations 
except that the field A^ plays the role of a multiplier field in the sigma model case (and 
hence has zero variation under Q). The twisted action can be written as the Q- variation of 
some gauge fermion A. The most general gauge fermion which is linear in the fermion fields, 
contains at most one derivative, and is not a Q-singlet is 



A = J d 2 xrj fl [ciA^ + c 2 d M + c 3 d u B^ lu ] 



(4.2) 



Furthermore, by a simple rescaling of the fields I can set c\ = \ an d c 2 = C3 = 1. Carrying 
out the Q-variation and integrating out the multiplier field A^ (along the imaginary axis) 
leads to the action 



S 



d 2 x 



As in [^5| the twisted fermionic action is of Kalher-Dirac form. Writing u 1 
bosonic action can now be rewritten 



11 



(4.3) 
2 - £10 the 



S B 



d 2 x 



1 

2 L 



2 



(4.4) 



where the implied summations over Roman indices run from i = 1,2. This action can be 
simplified by introducing the projection operator 



P +iu = - (5 i 5 v + e i e v ) 



(4.5) 



where the distinction between upper and lower indices is, for the present, immaterial. The 
bosonic action now reads 



S B 



d 2 x 



P^dvu 3 



(4.6) 
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Proceeding in this way we can introduce a new grassman valued field which is the superpartner 
to u % by defining ip 1 = A and = Xi2- Similarly the original field A„ (and its partner 77^) 
can be promoted to a field indexed by a Roman target index —* A~^ 1 provided the field is 
required to be self-dual under the action of this projector i.e 

4 = A* = P+?Ai (4.7) 

In terms of these new variables the Q-transformations are now 

Qu* = V 1 = , , 

Q V p = Af QAf = 1 ' j 

Going to complex coordinates we recover our original sigma model variations for a flat two 
dimensional target space. It should be clear that the sigma model fields are essentially the 
the Hodge self-dual components of the original Kahler-Dirac fields. It should also be clear 
how to proceed to a more general target space. Clearly one must be able to introduce a 
tensor Jj, here just e*-, which squares to minus the identity and is covariantly constant (so as 
not to disturb the form of the (^-transformations when applied to the self-dual fields). These 
requirements require the target manifold be Kahler with Jj its complex structure. Of course 
these were just the restrictions found in earlier constructions of the twisted sigma models [0]. 
Keeping the simple form of the (^-transformations the twisted sigma model action can be 
obtained as 

m ( 4 -9) 



S = Q j d 2 xrf dpu? - \a+i - 



Having recast the continuum theory in this language we can return the problem of con- 
structing a corresponding lattice model. We see an immediate problem. To avoid spectrum 
doubling we must discretize continuum derivatives with care. Specifically an exterior deriva- 
tive must be replaced by a forward difference operator D + in the discrete theory while the 
adjoint of an exterior derivative by a backward difference operator D~ if the resulting lattice 



theory is to be free of spectrum doubling [27]. Thus the bosonic part of the action in eqn. (| 
must be replaced with 

Sb = Y,\( D+ ^ + D ~ B ^) 2 
Since the difference operators are no longer the same we can no longer introduce the projector 
P + and rewrite the theory in terms of self-dual fields. Thus the Kahler-Dirac approach 
cannot be used to construct twisted lattice theories without spectrum doubling. A similar 
phenomena is encountered in the Yang-Mills case where the self-dual nature of the fields in 
four dimensional N = 2 super Yang-Mills prevents construction of a Q-exact and double free 
lattice theory. 

5. Simulations 

Here we revisit the 0(3) sigma model discussed in detail in (T|. The 0(3) sigma model has 
the advantage of being simple and can be reformulated as a twisted model with a CP 1 target 
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space. The metric, connection and curvature take the form, 



9uu = t 5 - 1 ) 



1 

2p 

in 

Ku = 9 uu du9uu = (5.2) 

P 

RuUUU 9uu9u^uU A (5-3) 

where p = 1 + uu. From ( |3.1| ), ( |3.4| ) and ([T^) the (effective) lattice action for a flat base 
space is given by (here 1 = 1 = 1), 



-2 uA 5 u + ^ (m W /u)(m W /n) - — -^-D^V ~ ^??£>+V' - ^jW^ 



p 2 + " — p 2- rr- /V "rr-, 2p2 ' -r 2p2 < +r ^ 

— 1 1 u u 



+ ln[/3/(2p 2 )] (5.4) 



z 1 u u 

~4^^2p^ mw + mH/ 27 " ^( m ^ n ) - 3(mKT!i)]i) 

where a factor of two has been absorbed into the coupling (5 and we have simplified our 

notation by replacing 1 — > u,ip l — > if), m_ — > and 77I — > 77. The explicit form of the lattice 

covariant derivative is _ 

£>£ = A* - — (A*u) (5.5) 
p 

Note that the new term ln[,5/(2p 2 )] missing in |TJ , emerges from integrating out the auxiliary 
fields B l + 2 . To proceed further it is convenient to introduce an other auxiliary field a to 
remove the quartic fermion term. Explicitly we employ the identity 

(3 V j Dae' a ^ + ^ + ^) = (5.6) 

where V is the number of lattice sites. Thus the partition function of the lattice model can 
be cast in the form 

Z = j DuDoDr)Di)e- s ^ a ^ (5.7) 
where the action is now given by 



1 . a . c 1 , w ,1 1 



^AluA^u+—(mwu)(m w u) + -era - -ln(2/> 2 ) + <SM(u,a)^ 
p z p z 2 p 



X 

where we have assembled the twisted fields into Dirac spinors 



(5.8) 



2 We thank Joel Giedt and Erich Poppitz for pointing out this factor 
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and the Dirac operator M(u, a) in this chiral basis is 



Mi, «\ - ( ^ mw ~ ^ mwu) + kc M A +-f ( A + n ) + 

M { u,a)- { {AS + ^ (A5 _ } _ _^ ^ mw _u {mwu)+hc 



(5.9) 



Notice the extra terms depending on the target space connection appearing along the diagonal 
which correspond to the twisted Wilson mass operator. These were not present in our previous 
paper Q. 

In order to be able to simulate this model, one needs to rewrite the effective action in 
a form that doesn't involve the grassman fields. For that we integrate out the Dirac field 
generating, 

det^ (/3 2 M(u, a)^M(u,a)) (5.10) 
The effective action is now given by, 

S = pS B (u,a) -ln(2p 2 ) - ^Tr In (/3 2 M f (u, a)M(u, <r)) (5.11) 

Clearly, the form of the fermion effective action we employ does not take into account any 
non-trivial phase associated with the fermion determinant - our simulation generates the 
phase quenched ensemble. We later examine the phase explicitly. 

To simulate this model we used the RHMC algorithm developed by Clark and Kennedy 
[ j3l[ l . The first step of this algorithm replaces the fermion determinant by an integration over 
auxiliary commuting pseudofermion fields F and F^ in the following way, 

deti(/3 2 M(n,(7)tM(n,a)) = J DFDF^"^ ^M)-h F (g 12) 

The key idea of RHMC is to use an optimal (in the minimax sense) rational approximation 
to this inverse fractional power. 

1 P{x) 

— r^j 

X 2 Q{x) 

where 

N-l N-l 

^ kx i 



1=0 8=0 



Notice that we restrict ourselves to equal order polynomials in numerator and denominator. In 
practice it is important to use a partial fraction representation of this rational approximation, 



i N 

+ (5.13) 

i=l 



- ' ' - 7~~T x + Pi 



The coefficient oti, Pi for i = 1,...,N can be computed offline using the Remez algorithm. 
Furthermore, the coefficients can be shown to be real positive. Thus the linear systems 
are well behaved and, unlike the case of polynomial approximation, the rational fraction 
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approximations are robust, stable and converge rapidly with N. The resulting pseudofermion 
action becomes 

' Ftp I VFt Ui 

' ' ^ M^M + f3i 



S PF p 



(5.14) 



In principle this pseudofermion action can now be used in a conventional HMC algorithm to 
yield an exact simulation of the original effective action |3^]. This algorithm requires that 
we compute the pseudofermion forces. For example, the additional force on the scalar field u 
due to the pseudofermions takes the form, 

i 

where the vector x % is the solution of the linear problem 

(M+M + Pi)x l = F 

The final trick needed to render this approach feasible is to utilize a multi-mass solver to 
solve all N sparse linear systems simultaneously and with a computational cost determined 
primarily by the smallest We use a multi-mass version of the usual Conjugate Gradient 



(CG) algorithm [33|. Note that for our simulations we have used N = 12 and an approxima- 
tion that gives us an absolute bound on the relative error of 8 x 1(T 6 for eigenvalues of M f M 
ranging from 10~ 6 to 10.0 which conservatively covers the range needed. 

5.1 Spectrum 

As a check of super symmetry, we have studied the boson and fermion propagators projected 
to zero spatial momentum, namely, G B (t) = Re(< u(t)u(0) >) and GfAt) =< \Pj (*)%(()) >. 
One can then extract the lowest lying mass states by fitting these two point functions to the 
form a + bcosh(m B (t — T/2)) for bosons and A(t — T/2)5-ij + iB(t — T/2)eij for fermions 
respectively. The quantities A and B are even and odd functions of their argument and 
we have taken the Dirac gamma matrix in the time direction to correspond to the usual 
Pauli matrix ai- To calculate the fermion masses rriQ , Wq, and mf we have used 
the simple functions acosh(m^ (t — T/2)), acosh(mf 1 (t — T/2)), asinh(m^ 1 (t — T/2)) and 
asinh(mf (t-T/2)) to fit ReG^ (t), ReGf^t), ImG^(t) and ImGf (t) respectively. Figure g 
shows the bosonic and fermionic masses versus coupling (3 in the range 0.5 to 10.0 for different 
lattice sizes. While the masses are quite different at small couplings, the mass degeneracy 
required by supersymmetry is recovered as we approach the continuum limit corresponding 
to large coupling (3. 

5.2 Ward Identities 

Of prime interest and a stringent test of supersymmetry are the supersymmetric Ward iden- 
tities. Consider first the scalar supersymmetry Q. The Ward identities corresponding to Q 
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are simply expectation values of the form < QO >. The simplest of these corresponds to the 
action itself, indeed as the latter is Q-exact Q, it is clear that, 



= ^ J [D^]e^ s QA =< QA > 



(5.15) 



which should be zero if the action is supersymmetric . On the other hand using the effective 



action in fl5,ll| ), 

- a -w=*=«*>-ir <-> 

Figure ^ shows a plot of < Sb > for a range of couplings (3 and different lattice sizes. 
While there are deviations from unity at small coupling these seem to disappear as the 
coupling (3 is increased beyond j3 ~ 4.0 and yield very strong support to exact Q-symmetry 
in the continuum limit. 

We have also looked at other Ward identities. Consider the local operator O of the form 

O = h + -g ll{x) W + {x)d^ 7 {y) + ijl{x)d+<f{y)] (5.17) 

The operator O is chosen in such a way that QO is local and leads to a 2-point function for the 
fermions. It is also Lorentz invariant in the base space and invariant under reparametrizations 
of the target space (at least in the continuum limit where the base difference operator becomes 
a true derivative). These constraints ensure the resultant supersymmetric Ward identity is 
satisfied non-trivially. In the 0(3) case < QO >= leads to 



< 



2p 



(5.18) 



Figures ||, |4] and |H| show plots of ( |5,18| ) projected to zero spatial momentum for different 
lattice sizes and coupling /? = 0.5,4 and 10.0 respectively. While for small coupling (3 the 
Ward identities are violated they are clearly satisfied for large coupling and vanishing lattice 
spacing. Notice the additional, at first sight rather startling feature; for (3 > 4.0 the functions 
are approximately independent of the coordinate separation \x — y\\ Actually this result 
follows from the Q-exactness of the theory. A continuum theory which is Q-exact has an 
energy momentum tensor T^ u = j^— which is also Q-exact. Furthermore, it then follows 
that the correlation function of two Q-invariant operators is independent of the metric and 



3 Note that, this tells us that Z is independent of /3, as long as /3 is not zero, and thus one can evaluate Z 
in the large-/? limit. Such a limit corresponds to the semi-classical approximation. Such an approximation is 
exact in our case (this is true for Witten type theories in general) 
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hence also of their separation |}4j. Using the Q- variation of the 77+ field in eqn. ( [2. 3D and the 
equation of motion B + = d + u we see that 

2u 

d + u = Qrj i/jt] (5.19) 

P 

The second term becomes small for large coupling (5 since u ~ f5~ 2 and hence in this limit the 
correlator is dominated by the leading Q-exact piece yielding the constant correlator as we 
observe. Thus the Ward identity reveals not only that the scalar supersymmetry is restored 
for large (3 but that the energy-momentum tensor and action of the theory are also Q-exact 
as expected from the continuum theory. 

We have also looked at the Ward identities following from the other supercharges of the 
twisted N = 2 supersymmetry (see Appendix [B] for the action of these other supersymme- 
tries). Indeed for the vector supercharges G± we have studied 

< h+~ (G+0_ + G-0+) >= (5.20) 



where, 



0_ = g ll {x)i, I {x)d.4> J {y) (5.21) 
0+ = g I j(x)ip 7 (x)d + ^ I (y) (5.22) 



For 0(3), eqn. (EM) gives 



1 _ 11 — 

< ^2(^ d + u ^ d - u ^ >= "2 < 2pHxj^ d+V ^ > (5 " 23) 



Figures ||, |7| and || show plots of eqn. 5.23j projected to zero spatial momentum for different 



lattice sizes and coupling. As for the scalar supercharge Q the Ward identities due to the 
vector supercharges are broken for small coupling but appear to be restored without fine 
tuning for small lattice spacing and large coupling. Again for large enough coupling these 
correlation functions appear independent of the temporal coordinate. This again follows from 
the fact that asymptotically we are examining the correlator of two Q-exact operators; Qrj 
and Qu = ip. 

We have additionally examined a simpler Ward identity. Namely 

\{ri{x)${y)) = {d + u{x)u(y)) (5.24) 

This Ward identity arises as a result of applying the scalar supercharge Q to the operator 
O' = r](x)(j)(y), i.e. {QO') = 0. Figure ^ shows a numerical calculation of eqn. ( |5.24| ) 
projected to zero spatial momentum. While supersymmetry is broken at small (3 coupling it 
is clearly restored as we approach the continuum limit, confirming the results we obtained 
with the other Ward identities. 
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5.3 Fermionic Condensate 



In the previous section we considered Ward identities that lead to fermionic 2-point functions. 
Here we consider Ward identities involving 4-point functions of the fermion fields. As we 
will see, the physics here is richer and more complicated as these quantities can receive 
contributions from instantons. Calculations in the continuum theory ^] predict the 
presence of a vacuum condensate < ^ip^ > similar to the gluino condensate in J\f = 1 super 
Yang- Mills theory. Since ifi = ^ l a single chiral component of the Dirac field we see that the 
presence of such a condensate constitutes an anomalous breaking of the chiral symmetry of 
the theory. 

To see how this comes about in the twisted theory consider the following correlation 
function 

C(x - y) =< G{x)Q(y) > (5.25) 

where 

9(x) = J j7 ^V(aO ( 5 - 26 ) 

Here, <// is the complex structure in the target manifold which locally takes the form iSj . 
It is straightforward to verify that 0(x) is invariant under Q. Hence, as we argued earlier 
the correlation function C(x — y) should actually be independent of \x — y\; C(x — y) = C 
a constant. Further, using cluster decomposition, it can be shown that any non-vanishing 
value for C implies a non-vanishing value for the condensate (Q(x)} = (0) = i-v/C*. 4 At first 
glance the classical chiral symmetry of the theory would prohibit such a non-zero value for C. 
However quantum anomalies can and do spoil this symmetry. As is well known, the theory 
admits non-trivial classical solutions called instantons. These are given by solutions of the 
equations 

d + u = (5.27) 
The simplest single instanton is hence given by the analytic function 

•M = 7^3 (5-28) 

where a and (5 are complex constants. There are four real zero modes associated to variation 
of these parameters and hence, by supersymmetry, there will be four real fermion zero modes 
localized in the vicinity of such an instanton. Examination of the equation D + \E' = in 
such a background reveals these zero modes to be chiral and such a condensation will hence 
break chiral symmetry. Furthermore, since is topological its expectation value can be 
computed exactly in the semi-classical limit corresponding to a one loop computation in such 
an instanton background. Notice that the number of fermion zero modes is just sufficient 
to saturate the four point function we are considering. The result is given in [35]. A lattice 



4 A word of caution here - (37], |3^] shows that cluster decomposition fails in the strong coupling instanton 
calculation of the analogous condensate in JV = 1,D = 4 super Yang-Mills theory. Here, the topological 
character of C allows it to be calculated exactly at weak coupling 
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calculation of this quantity is potentially very interesting as it yields information on both 
the U(l) anomaly, the presence of approximate lattice instantons and, as we will see, the 
dynamical breaking of super symmetry. 

In our lattice regulated 0(3) model we have hence measured the four point function 

< >=< > ( 5 - 29 ) 

P (°) P w 

The plots in figure 1C. shows ( |5,29| ) projected to zero spatial momentum for different lat- 
tice sizes and different values of the coupling (3. Clearly in the large f3 limit < @(t) > is 
approximately independent of t confirming the presence of a condensate. 

Of course on a torus only instanton /anti-instanton pairs can exist but provided these are 
well separated and the discretization errors small we can still expect a condensation of four 
approximate zero modes in the vicinity of such a lattice instanton yielding a corresponding 
contribution to the local value of C. In practice the value of the correlator is not constant as 
we will find contributions to C also from the anti-instanton. This is particularly true for small 



lattice volumes as can be seen in the figures 10. There are also explicit lattice supersymmetry 



breaking effects visible in the data at f3 = 0.5. Notice that the value of the condensate for a 
fixed lattice volume initially grows as (3 is increased but eventually turns over and decreases. 
We can understand this as the result of finite volume effects which become large for large (5 
and suppress such instanton-like configurations. Notice that the correlation function we plot 
is dimensionless and hence contains powers of the lattice spacing a((3). We have not measured 
this quantity directly so it hard to assess from figure ^ whether this condensate survives the 
continuum limit. We will try to address this question in the next section. 

We have argued that the presence of such a condensate is signal for the anomalous break- 
ing of chiral symmetry but it can also be viewed as an order parameter for supersymmetry 
breaking. To see this notice that the condensate of the 0(3) model can actually be obtained 
as the Q-variation of the following operator. 

Ocond = , ; 2 ; -M x)i>(y)xlj(y) (5.30) 

Thus a non-zero value for this quantity is equivalent to the statement < QO >^ and 
indicates a dynamical breaking of supersymmetry driven by instantons rather along the lines 
originally envisaged by Witten [39]. However it was pointed out in [35, [3(| that this operator 



Ocond is not invariant under the global 0(3) symmetry of the model. It is conjectured that 
supersymmetry is only violated in the unphysical 0(3) non-invariant sector of the theory. 
Specifically, the expectation is that any operator QO containing a multiple of a four-fermion 
term and for which O is not 0(3) invariant would develop a vacuum expectation value different 
from zero whereas for an 0(3) invariant O it would remain zero. For a general N = 2 sigma 
model the notion of 0(3) invariance would then be replaced by the requirement that the 
operator be invariant under all possible isometries of the target metric. It would be interesting 
to investigate these issues in more detail. 
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5.4 Phase 



Up to this point we have neglected a possible phase arising in the calculation of fermion 
determinant. It is not hard to show that the determinant is real since generically the eigen- 
values occur in complex conjugate pairs. The exception occurs with purely real eigenvalues. 



As one of these crosses the origin the sign of the determinant will change. Indeed figure [11. 
indicates that this phase does undergo weak fluctuations and one has to check for its effect 
on the simulation. As usual we can always compensate for neglecting the phase factor in the 
simulation by re-weighting all observables by the phase factor according to the simple rule 

< o >= (5-31) 

Hence we examined the Ward identities we have studied in the previous section now weighted 
with this phase factor. For example table |l] shows the mean re-weighted bosonic action 
together with the phase quenched value for different lattice sizes and a range of couplings 
(3. While re-weighting typically amplifies the estimated errors it does not appear to change 
the mean value for this observable at least within statistical errors. This seems to indicate 
that the phase fluctuates approximately independently of the remaining part of the measured 
observable leading to an, at least approximate factorization in the reweighted observable 

_^ < Oe ia > Q=0 < Q > a=0 < e ia > a=0 

< O > f uii= ■ . =< O > a=0 (5.32) 

< e ia > a=0 < e ia > Q=0 



5.5 Continuum Limit 

Figure [l2| shows a plot for the dimensionless ratio value <1// ^ t ^ > (m b is the lightest boson 
mass) for different lattice sizes and for different values of the coupling. One can see that 
an approximate plateau occurs for sufficiently small lattice spacing which we interpret as 
the onset of continuum physics. The separation of the curves from different lattice volumes 
we consider to be a measure of possible finite volume effects. Modulo these finite volume 
questions we take this as evidence that the magnitude of the condensate is non-zero in the 
continuum limit. 



6. Conclusions 

This paper is devoted to a study of the N = 2 supersymmetric sigma model regulated on 
a lattice. The lattice formulation we employ was derived in an earlier paper and results 
from a discretization of a twisted version of the continuum theory. The twisting process is 
to be viewed as a change of variables (we are in flat space) and has the merit of exposing a 
nilpotent scalar supersymmetry Q which can be preserved under discretization. In this paper 
we extend our previous work by introducing a modified Wilson operator in the form of a 
twisted mass term which is used to remove doubler modes from our lattice theory without 
spoiling the Q-exact property of the lattice action. The introduction of such a mass term is 
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p 


P<s b > 


< Sb >W 


0.5 


145.523 ± 0.286 


144.464 ± 3.858 


1.0 


147.369 ± 0.281 


148.764 ± 21.547 


2.0 


135.51 ±0.281 


139.749 ± 53.707 


3.0 


128.275 ± 0.242 


122.731 ±8.747 


4.0 


127.272 ± 0.329 


126.933 ±2.093 


5.0 


127.326 ± 0.213 


127.342 ± 1.274 


10.0 


127.552 ± 0.299 


127.164 ± 1.967 



P 


(3<S b > 


(3 < Sb >w 


0.5 


321.541 ±0.568 


321.461 ± 8.984 


1.0 


310.143 ±0.391 


305.595 ± 17.114 


2.0 


291.275 ± 0.360 


284.338 ± 40.365 


3.0 


286.719 ± 0.352 


276.315 ± 24.448 


4.0 


285.546 ± 0.326 


283.904 ± 5.489 


5.0 


285.627 ± 0.340 


285.307 ± 3.386 


10.0 


287.163 ± 0.285 


286.943 ± 2.483 



(3 


P< s b > 


(3 < Sf, >w 


0.5 


283.02 ±2.166 


285.369 ± 4.369 


1.0 


458.823 ± 1.707 


459.028 ± 7.847 


2.0 


509.266 ± 0.628 


504.264 ± 38.233 


3.0 


503.627 ± 0.784 


496.215 ± 46.938 


4.0 


507.243 ± 0.609 


505.991 ± 13.837 


5.0 


507.412 ± 1.028 


506.915 ± 14.778 


10.0 


511.561 ±0.570 


511.715 ±5.386 



Table 1: (3 < Sb > for quenched and unquenched phase on 8x8, 12x12 and 16x16 lattices. 



compatible with supersymmetry if it takes the form of a holomorphic Killing vector in the 
target space and corresponds to the addition of central charges in the supersymmetry 
algebra. The replacement of a simple mass term with a Wilson difference operator breaks 
supersymmetry softly and is not expected to lead to additional fine tuning. We study the 
CP 1 ~ 0(3)-model in detail using the RHMC algorithm on lattices as large as 16 over a 
range of coupling (3 = 0.5 — ► 10.0. Our results for both the spectrum and Ward identities 
provide strong evidence for a restoration of full supersymmetry at large f3 and small lattice 
spacing without additional fine tuning. Notice that while the Ward identities are examples 
of trivial topological observables the mass spectrum we measure is non-top ological. Further 



investigations of such physical observables are underway and will published elsewhere [41]. 
We are also able to see a non-zero chiral condensate and this appears to persist into the 
continuum limit as expected from the chiral anomaly. It would be very interesting to extend 
this work to the general CP N ~ l models where the numerical results could be compared with 



exact results on the low lying mass spectrum [42]. 
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A. Appendix: Notations 

A real system is denoted by x , x 2 which combine to give holomorphic coordinates, 

z = x +ix 2 , z = x l — ix 2 (A.l) 
Then a vector with components (V 1 , V 2 ) has holomorphic components, 

V + = hp x - iV 2 ), V. = ~(V! + iV 2 ) (A.2) 

The Dirac matrices 7^ are given by (Chiral basis), 

(7V = ^\ (7V = ^ 2 (A.3) 
Spinor indices are raised and lowered by the matrix C a p = a , 

W)afi = (injCrfi (A.4) 

thus for a vector V^, 

it+V^ = Vi- iV 2 = V + (A.5) 

1 »_V fl = Vi + iV 2 = V- (A.6) 

(A.7) 

For a Kahler manifold, the metric is given by, 

9lJ - = djdjK (A.8) 

where K is a Kahler potential. (Note that gjj = gu = for a almost complex manifold). 
The non vanishing components of the Christoffel symbols are then, 

^jk = 9 lL dj9KL, ^jk = 9 lL dj9 L K (A.9) 
This implies the only non-trivial components of the Riemann tensor are, 

Rijkl = 9im9k^jl, or R 1 jRl = d^T^, (A.10) 

B. Appendix: Twisted M = 2 supersymmetry algebra 

The algebra of N = 2 supersymmetry in 2 space dimensions is given by [28], 



{Qa+, Qp-} = laP P ^ i R ' Q<*±] = ±^Qa± 

{Qa+,Qp+} = {Q a -,QfS-} = 0, [J,P„] = -ie»P u 

[Qaa, Pp] = [Pp, PA = 0, [R, P^] = 

[J, Q± a ] =±\Q ±a , [J, R] = [J, J] = [R, R] = 



(B.l) 
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where J is the generator of Lorentz 50(2) symmetry and R is the generator of the internal 
SO (2) symmetry. In order to obtain some scaler charges, one perform a twisting that changes 
the spin of the above charges. In order to do this we need to redefine the Lorentz generator. 
Indeed if we take (see Appendix [A| for our notations), 

J = J + R (B.2) 

with respect to this new generator one finds that and Q (_ behave as scalars while the 

pair Q++ and Q as vectors. To make manifest the new Lorentz structure of each of the 

generators we define, 

Ql = Q+- 
Qr = Q-+ 
Q ++ = ^ + G, = G + 
= j»_G„ = G_ 



it is clear from flB.l| ) that, 

Ql = Qr = {Ql,Qr} = (B.3) 



The algebra ( |B.1| ) in terms of the new Lorentz generator J and the following redefined oper- 
ators, 

Q = Ql + Qr, M = Q l -Q r (B.4) 

is 

{Q,G fl } = P ll 

[Q, J] = [M, J] = 

{M,G^} = -ie l /P v (B.5) 

[J,P„} = -ie f ?P I/ 

[J,G,] = -ie"G v 

[P^P V ] = {G^G V } = [J, J] =0 

in addition, the action of the R generator on the new twisted charges is, 

[R,Q] = -\M 

\RM\ = -\Q v (B6) 

[R,Gfj\ — ~2 e [l G u 

[R,R] = [R,J] = [R,P„] = 

Let our field content be 4> l , B l a ,ip l and rf a or on the complex manifold cj) 1 , B^i/j 1 and rjL ■ 
The R-transformations of these fields is given by 

/Ji _ n r D jJl 



[R,P]=0, [R,P]=0 

[R,^ I ] = ±i> 1 , [R,^} = -^ 

[R,v{] = W+> [R,vi)=-h- ' 

[R,B I + ]=B I + , [R,BL]=-BL 
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More importantly the transformation of our fields under the twisted charges Q,M, G+ 
and G_ is given by, for Q, 

Qcp 1 = i) 1 Qcf) 1 = V J 

Qtf = QV 7 = 

Qvi = B{~ T kj^ J V+ QV 1 - = Bl- r^ J ??5 1 ' ' 

QBl = -Y^^Bf - R^j^^ri QBi = -T ! m iP J B* - B 7 jLR i, J ^ L r, R 

for M 

Mcj) 1 = -ip 1 Mcf) 1 = i\) J 

Mtf = MV J = 

Mr}\ = -{B\- V KJ ^ J r}^) MrjL = Bi - r^ J r?5 ' 

MB J + = Tj K ip J B^ - R 1 Klj ^ J ^ L rft MBl = -T JR 7p J B R + R 1 RLj ip J ip L ri* 

for G+ 

G+cf) 1 = \r)\ G + <f> r = 

G + r] J + = G+T]L = 

G+Bl = -iT^BinX G+Bl = d +V L - Wklj^^- 
and finally for G_, 

= G-4> 1 = \rjL 

G-ip T = G_^ f = -\{B ] _ -T R3 ^ J r] R ) 

G-r ] I + = G-7]l=0 
G-B{ = d-i]{ - \R : K - LJ ^ J v L rif G-Bi = -\T J 3R Blrf_ 



(B.10) 



(B.11) 
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Figure 1: mass spectrum for 8x8, 12x12 and 16x16 lattices 
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Figure 2: Simple Ward Identity 
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<QO>=0 

8x8 lattice, p=0.5, O^g^fxJl^d.ubarj^bar^uy 



> bosonic part 

> fermionic part 




<QO>=0 

8x8 lattice, p=4, 0=g uubar (x)[i|j > (d_ubar) y +x|/bar > (d i u) y ] 



> bosonic part 

> fermionic part 




<QO>=0 

8x8 lattice, p=10, O=g„ 0b Jx)[\|;,,(d_ubar) y .H|/bar,,(d,uy 



bosonic part 
fermionic part 




Figure 3: < QO >= for 8 x 8 lattice 
- 25 - 



0.02 



<QO>=0 

12x12 lattice, p=0.5, 0=g uubar (x)[ v>; (cl_ubar) y + > );bar x (d t u),] 



bosonic part 
fermionic part 



0.01 - 



-0.01 



0.01 



0.005 




15 



<QO>=0 
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Figure 4: < QO >= for 12 x 12 lattice 
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Figure 5: < QO >= for 16 x 16 lattice 
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Figure 6: < h+-(G+0- + G-0+) >= for 8 x 8 lattice 

- 28 - 



<h + "(G + O_+G_O>=0 



12x12 lattice, (3=0.5, O t =g uubar (x) v>: (d_ubar) y , 0_=g uubar (x) v bar(d t u) ( 
0.01 

-o bosonic part 
-• fermionic part 



0.005 



-0.005 



-0.01 



10 



15 



<h + "(G + O_+G_O + )>=0 

12x12 lattice, (3=4.0, O t =g uubar (x) Vx (d_ubar) ( , 0_=g uubar (x) v bar(d t u) ( 



0.01 



0.005 - 



-0.005 - 



bosonic part 
• fermionic part 



-0.01 



10 



15 



<h + -(G + O_+G_O>=0 



12x12 lattice, |3=10.0, O^g^JxJv^d ubar) , O =g uuba ,(x)ybar(d t u) y 
0.005 

> bosonic part 
• fermionic part 



0.0025 - 



-0.0025 




Figure 7: < h + ~(G + 0- + G-0+) >= for 12 x 12 lattice 
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Figure 8: < h+-{G+0- + G-0+) >= for 16 x 16 lattice 
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Figure 9: < Q(r/ X u y ) >= for 8 x 8 lattice 
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Figure 10: condensate for 8x8, 12x12 and 16x16 lattices 
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Figure 11: Distribution of phase on 16x16 lattice 
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